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Abstract 

We propose a general and compact scheme for the computation of the periods and am- 
pUtudes of the chiral persistent currents, magnetizations and magnetic susceptibilities in 
mesoscopic fractional quantum Hall disk samples threaded by Aharonov-Bohm magnetic 
field. This universal approach uses the effective conformal field theory for the edge states in 
the quantum Hall effect to derive explicit formulas for the corresponding partition functions 
in presence of flux. We point out the crucial role of a special invariance condition for the 
partition function, following from the Bloch-Byers-Yang theorem, which represents the 
Laughlin spectral flow. As an example we apply this procedure to the parafermion Hall 
states and show that they have universal non-Fermi liquid behavior without anomalous os- 
cillations. For the analysis of the high-temperature asymptotics of the persistent currents in 
the parafermion states we derive the modular S -matrices constructed from the S matrices 
for the m(1) sector and that for the neutral parafermion sector which is realized as a diagonal 
affine coset. 
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1 Introduction 



Mesoscopic physics has recently become a very intensive area of research because 
of the exciting opportunities it offers to observe fundamental implications of quan- 
tum mechanics in experimentally realizable two-dimensional samples. One impor- 
tant characteristics of the mesoscopic phenomena is that the electrons are coherent 
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which implies that the sizes should be extremely small, the samples extremely pure 
and with high mobility and the temperatures very low. 



According to a famous Bloch theorem the free energy of a conducting ring threaded 
by magnetic field is a periodic function of the magnetic flux through the ring with 
period one flux quantum (po = h/e = I. The flux dependence of the free energy 



where Z(T, (p) is the partition function gives rise to an equilibrium current 



called a persistent current, flowing along the ring without dissipation. These oscil- 
lating currents have universal amplitudes and temperature dependences and have 
been observed in mesoscopic rings [1], where the length of the ring is smaller than 
the coherence length at very low temperatures. 

The fractional quantum Hall (FQH) states could be genetically interpreted as meso- 
scopic rings because the only delocalized states live on the edge of the two-dimen- 
sional disk sample. The universality classes of the FQH systems have been success- 
fully described by the effective field theories for the edge states in the thermody- 
namic scaling limit [2]. These turned out to be topological field theories in 2 -i- 1 
dimensional space-time which are equivalent to conformal field theories (CFT) on 
the 1 -I- 1 dimensional space border, provided the bulk excitations are suppressed by 
a nonzero energy gap. 

Most of the computations of the persistent currents in strongly correlated meso- 
scopic systems, including in the FQH regime, so far have been based mainly on the 
Fermi or Luttinger liquid picture [3,4] because a more general approach to these 
quantities was missing. In this paper we shall formulate a general and compact 
scheme for the computation of the chiral persistent currents in mesoscopic FQH 
disk states based on the rational CFT on the edge. This scheme has become very 
important after recent experiments [5] showed that the (multicomponent) Luttinger 
liquid theory is irrelevant for more general filling factors such as those in the prin- 
ciple Jain series. We stress that the rational CFT description which is the right strat- 
egy for classifying the FQH universality classes is also the right starting point for 
the computation of the persistent currents in the FQH states. The main advantage 
of this approach is that the CFT partition function which is used as a thermody- 
namical potential can be computed analytically even when Aharonov-Bohm (AB) 
magnetic field is threading the sample. This fact has already been used in Refs. [6] 
and [7] for the computation of the persistent currents, magnetizations and the cor- 
responding magnetic susceptibilities with the intuitive explanation that because the 
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transformation ^ ^ ^ + r represents the addition of one flux quantum its general- 
ization ^ ^ ^ + (f>T should be interpreted as adding flux (p in units (f>o = I. We derive 
this result in Sect. 2.8 which was one of the motivations for the present paper. The 
fact that the magnetic field is of AB type (i.e., an infinitely long and thin solenoid) 
is a significant simplification because the AB magnetic field is zero at the edge so 
that it does not break the conformal invariance. On the other hand this limitation is 
not completely restrictive because the edge of a FQH system is so thin that the free 
energy depends only on the flux encircled by the edge and not on the magnetic field 
alone. Adding AB flux to the FQH system naturally introduces twisted boundary 
conditions for the field operators of the electron and quasiparticles and also modi- 
fies the efl'ective Hamiltonian. We show that the ultimate eff'ect of this twisting on 
the partition function is that it is multiplied by a universal exponential factor and 
that the modular parameter f is shifted by an amount proportional to the modular 
parameter t where the coefficient of proportionality coincides with the AB flux 0. 
As an illustration of the general CFT approach to the chiral persistent currents we 
are computing these currents for the parafermion FQH states using the partition 
functions developed in [8]. 

The persistent currents in the FQH regime typically contain large non-mesoscopic 
components [9], slowly varying with the magnetic field, which cannot be described 
within the CFT framework since there are contributions from states laying deep 
below the Fermi level. The edge states effective CFT is relevant only for the com- 
putation of the much smaller oscillating part of the persistent currents which could 
nevertheless be tested directly in SQUID detectors [1] using the Josephson effect. 
While the persistent currents naturally express the electric properties of the sys- 
tem, it turns out that they could also give important information about the neutral 
structure of the effective field theory for the edge states [3,10]. As we shall demon- 
strate in Sect. 4.3 and Sect. 4.4 the low- and high- temperature decays of the per- 
sistent current's amplitudes crucially depend on the neutral properties of the edge 
states. Our analysis reveals two different universal mechanisms for the thermal re- 
duction of the persistent currents: at low temperatures the thermal activation of 
quasiparticle-quasihole pairs leads to a logarithmic decay of the amplitudes while 
at high temperatures the reduction is due to thermal decoherence and the decay is 
exponential. 

The rest of this paper is organized as follows: in Sect. 2 we review the rational CFT 
approach to the FQH states, describe in brief the terminology of the chiral quan- 
tum Hall lattices of topological charges which is necessary for the understanding 
of the general classification of FQH universality classes. In Sect. 2.8 we show how 
the CFT partition function is modified by the AB flux starting from the twisting of 
the electron field and deriving explicitly the modification of the m(1) current and 
the Hamiltonian. In Sect. 2.9 we discuss the Laughlin spectral flow, the role of the 
Cappelli-Zemba factors for restoring the invariance of the partition function un- 
der this flow and prove that the CFT partition function with AB flux is invariant 
under ^ + 1 as it should be according to the Bloch-Byers-Yang theorem. In 
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Sect. 2.10 we derive the low-temperature asymptotics of the CFT partition func- 
tion, which is universal. This asymptotics is then used in Sect. 4.3 to derive the 
low-temperature asymptotics of the persistent currents amplitudes which have a 
universal form and is expected to be valid for all FQH states. In Sect. 2.1 1 we show 
how to use the modular S transformation in a rational CFT to analyze the high- 
temperature asymptotics of the partition functions and the persistent currents. In 
Sect. 3 we summarize the results from Ref. [8] about the partition functions for 
the edge states of the parafermion (PF) FQH states for k = 2,3 and 4. We de- 
rive analytically in Sect. 3.2 the complete S matrix for the parafermion states as 
well as that for the neutral part which is realized as a diagonal affine coset. Then 
we investigate the persistent currents for fixed temperature when the flux is varied 
within one period as well as the temperature decay of the amplitude. In Sect. 4.4 
the high-temperature asymptotics for the parafermion FQH states is computed us- 
ing the analytic S matrices derived in Sect. 3.2. Some technical details are put in 
several appendices. 



2 Chiral partition function for the FQH states with Aharonov-Bohm flux 

2.1 CFT for the edge states: labeling the FQH universality classes 

Due to the existence of a nonzero energy gap in the FQH effect excitations spectrum 
the electron system behaves as an incompressible electron liquid at low energies. 
The effective field theory in the scaling thermodynamic limit, which essentially de- 
scribes the low -energy excitations, turns out to be a topological relativistic quantum 
field theory in 2-1-1 dimensions with a Chem-Simons action. It is well known that 
these theories are equivalent to CFTs on the 1 + 1 dimensional edge. Therefore, it 
is believed that the CFT could capture all universal properties of the FQH systems 
and thus allow us to make a complete classification of the QH universality classes. 

The electric current on the edge, where the coordinates are denoted by {Xa} for 
a = 1, 2 and {A^} is the restriction of the vector potential to the edge, is anomalous 
[11] 

1 ^2 

i.e., it violates the m(1) invariance on the edge, however, this is exactly compensated 

by the violation of the m(1) invariance in the bulk so that the sum of the bulk and 
edge current is invariant. This compensation of bulk and edge anomalies is the ori- 
gin of the so called holographic principle [11] which in this case can be formulated 
as the correspondence between the quasiparticle spectra on the edge and in the bulk. 
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The m(1) anomaly mentioned above is the basic mechanism through which the Hall 
current is transmitted between the two edges of the real FQH samples. 

2.2 The m(1) X Vir symmetry 

The chiral CFT describing the edge of a FQH system always contains a separate 
m(1) and Virasoro {Vir) symmetries. The first one describes the electric charge Q, 
or equivalently the magnetic flux properties due to the fundamental FQH relation 
[2] 

!2 = vh0, (3) 

{e = \) while the second corresponds to the angular momentum. Note that due to 
the Sugawara formula the m(1) sector contributes to the angular momentum as well 
so by Vir in the relative u{\)xVir decomposition we mean the Virasoro contribution 
from the neutral sector only. As we shall see in Sect. 2.4 the neutral Vir sector must 
be present for all FQH states with filling factors whose numerators are bigger than 
1. 

The chiral partition function Z{t, corresponding to a disk FQH sample with a 
single edge computed within the eff'ective CFT for the edge states 

N 

Z(T, = tr ^^-^/24 (2ni^j^') = YjXAir, 0, (4) 



is the linear sum of all topologically distinct chiral CFT characters [7] 

XA(r,0= tr q"^-''^^ exp(2m^4) where 'H = ® 'H;^. (5) 

"H^ ^ ' /l=l 

The number N in Eqs. (4), called the topological order [12], (the number of topo- 
logically inequivalent quasiparticles with electric charge [11] < 2ei < 1) is one 
of the most important characteristics of the FQH universality classes. The Hilbert 
space of the edge states H, over which the traces in Eqs. (4) and (5) are taken, 
is the direct sum of all independent (i.e., irreducible, topologically inequivalent) 
representations 'Ha of the chiral algebra which are closed under fusion (see the 
Introduction of Ref. [13] for a short description of the chiral algebra terminology 
in the FQH efi^ect context). The efi'ective Hamiltonian for the edge states in the 
thermodynamic limit is given by [14] 

vf(^ c\ , £ dz ^ (vpt-ixX 
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where x is the coordinate on the edge, t is the imaginary time, vp is the Fermi 
velocity on the edge, R is its radius, Lq is the zero mode of the Virasoro stress 
tensor T(z) and c is the central charge of the latter [15]. The CFT Hamiltonian (6) 
captures the universal properties of realistic FQH systems and could eventually give 
a qualitative description of the Hilbert space structure of the low-lying excitations 
[16], the energy spectrum and even the (universal part of the) energy gap [13] of the 
FQH systems in the thermodynamic limit. In particular, it should not be surprising 
that the CFT-based computations [10] of the persistent currents in the Laughlin 
FQH states give the same results as those based on the Luttinger liquid picture [4]. 

The electric charge operator in Eq. (4), defined as the space integral of the charge 
density J^\z) on the edge, which is proportional to the normalized m(1) current 

J(z) = i d4>^'\z) = J]jn z-"-\ [J„, Jm\ = n 5„+^.o, (7) 

n€Z 



has the following short-distance operator product expansion [2,15,14,13] 

Jt = ip-. J'\z), r\z) = V^7(z), r\z)r\w) ~ (8) 

J 2m (z - wY 



The modular parameters q and ^ entering Eq. (5) are generically restricted in annu- 
lus samples for convergence by 

q = Q^", ImT>0, {\m.T -llkeT) ^gC 



and their real and imaginary parts are related to the inverse temperature, chemical 
potential and Hall voltage [14]. For the disk FQH sample, however, both r and ^ 
have to be purely imaginary in order for the partition function Eq. (4) be real. The 
physical meaning of Imr and Im^ will be discussed in Sect. 2.9. 

Because the AB flux modifies only the charged m(1) sector of the model it would 
be convenient to artificially split the stress tensor into a charged and a neutral parts 

nz) = T^'\z) + T^''\zX where T^'\z) = ]^:J{zf: (9) 



is the Sugawara stress tensor of m(1) and T^^^ = T - T^'^^ is its complement. 
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2.3 Abelian rational CFTs: FQH states classification in terms of integer lattices 
of topological charges 



The CFTs which naturally appear [2] in the FQH effect context are rational exten- 

-N 

sion of the u(l) current algebra in terms of integer lattices of topological charges. 
A central result in this approach is the complete A-D-E classification of the FQH 
universality classes in terms of abelian rational CFT associated with those lattices 
[2]. We recall that an abelian theory is a rational conformal field theory with in- 

N 

teger central charge c = N and m(1) affine symmetry; the current algebra is gen- 
erated by abelian currents = id^(p"(z), extended by the vertex operators 
Y(A,z) = :exp(iAa(p"(z)y., whose "topological charges" A,, form a A'^-dimensional 
lattice F. The lattices which appear in the quantum Hall effect must satisfy some 
specific physical conditions which are summarized bellow. Following Ref.[2], we 
call chiral quantum Hall lattice an odd integral lattice to which a special charge 
vector Q from the dual lattice F* is associated. This vector assigns to any lat- 
tice point q e Y the electric charge {q\Q) of the corresponding edge excitation. 
In particular, the lattice should contain an electron excitation q with unit charge 

= The norm \q\-^ = {q\q) of the vectors ^ € F gives twice the confor- 
mal dimensions, i.e. twice the spin of the corresponding excitation, while the scalar 
product {q\q') is the relative statistics of the excitations q and q' . 

The charge vector Q satisfies the following conditions: 

(i) it is primitive, i.e., not a multiple of any other vector q* G F*; 

(ii) it is related to vh by: 

\Q\^ = vh = '^, nH,dH^n, gcd{nH,dH)=h (10) 
— dH 

(iii) it obeys the charge-statistics relation for boson/fermion excitations: 

(Ql£) = l£p mod 2 , for any £ G F. (11) 

We shall briefly recall the notion of maximal symmetry [2] since it covers the largest 
class of QH lattices. A chiral Hall lattice is called maximally symmetric if its neutral 
sublattice (Fq ± Q) coincides with the (internal symmetry) Witt sublattice Y-^ <zT 
(defined as the sublattice generated by all vectors of square length 1 or 2) so that 
dim Fjy = dim F - 1. The Witt sublattice is always of the type A ® D ® E ® Z. 
The maximally symmetric lattices, denoted in [2] by the symbol (L| -Fv^^), can be 
characterized in an appropriate basis {^'} by the following Gram matrix G = (G'O 
where G'^ = (q'\q^) is the metrics of the lattice: 



G = 



L 


CO 


T 

or 


CiYw) 



det(G) = — det [CiTw)] ; (12) 
Vh 



1 



the charge vector is Q = (1, 0, . . . , 0) in the corresponding dual basis. Here L e N 
is the minimal relative angular momentum ^ that appears as the minimal power of 
{Zi - Zj) in the corresponding wave function; C(rTy) is the Cartan matrix for the 
Witt sublattice and a; is an admissible weight for restricted by the condition 
< L (see Eq. (5.6) in [2]). Note that the Witt sublattice does not contribute to 
the filling fraction vh since it is orthogonal to the charge vector Q. It is easy to find 
Vh corresponding to the lattice with the Gram matrix (12) using Eq. (10), 

Vh = (Q\Q) = Q^.Gr'.Q = , ] . , • (13) 
L - (ojIoj) 



The non-equivalent quasi-particles in the Hall state correspond to the irreducible 

representations of the extended algebra and are labelled by the points A e Y* /Y 
where A is the quasiparticle's topological charge and Y* is the lattice dual to Y. We 
recall that Y*/Y is a finite abelian group of order (detP) whose multiplication law 
represents the fusion rules. There are several important consequences of the specific 
structure of the Gram matrix (12) and FQH eff'ect specifics which we consider in 
more detail in the next subsections. 



2.4 The electron charge decomposition 



The Gram matrix (12) implies that the topological charge q of the electron should 

—el 

take the form [2] 

q=-—Q + ai, (14) 
^1 Vh — 



where Q is the lattice charge vector and co, which is the weight entering Eq. (12) 
i.e., an admissible weight of the Witt sublattice [2], represents the neutral degrees of 
freedom of the electron because {Q\co) = 0. Then the electron field operator would 

be simply a tensor product of a m(1) vertex exponent of a chiral boson 0, which is 
a simple current of the m(1) theory, and in general a simple current from the 
neutral part of the theory 

Mz) = :exp |-j-^^(^)(z)j: ® 0,(z), (15) 



where (p^'^'iz) is the normalized chiral m(1) boson corresponding to the charged com- 
ponent of the electron and f^wiz) is a primary field (oj is a simple current, i.e., 

^ In this case, it is equal to twice the conformal dimension of the electron operator, i.e. an 
odd integer. 
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co*co = l [15,17]) from the neutral sector with CFT dimension (a)\co)/2 
<P^'W\w) ~ - ln(z - w), ^M)K(w) ^ 



Next, due to the charge-statistics relation (11), the conformal dimension of the 
electron which is 

^«i = ^(£eilli> = ^(^ + (^i^^ = I^ + ^^"^(^i^' A^^o^ = (ojIo;) > (16) 



must be equal to p + 1/2 for some p e Z+. According to Eqs. (14) and (16) the 
neutral component (o of the electron should certainly give a nontrivial fractional 
contribution to the CFT dimension when the numerator of the filling factor vh is 
bigger than 1. Therefore, in general, the neutral component 0^(z) of the electron 
field operator cannot be completely removed without violating the fermionic statis- 
tics of the electron. Note that the condition (a)\oj) > in Eq. (16) is equivalent to 
/max > l/vff as formulated in Ref. [2]. 



2.5 Spin-charge separation and Pairing rule 



Because the weight coi^O has at least one non-zero component when the numerator 
of the filling factor is Hh > 1 the Gram matrix (12) in this case is not decompos- 
able, i.e, it cannot be represented as a direct sum of two matrices with smaller 
dimensions. This implies that the dual lattice would not be decomposable either 
so that the inequivalent irreducible representations of the chiral algebra, i.e., the 
quasiparticles, would not be simple tensor products of charged and neutral vertex 
exponents but actually sums of such. Therefore the charge and spin are entangled 
and this absence of spin-charge separation makes the situation more complicated. 
Nevertheless, the spin-charge decomposition could be achieved at the expenses of 
enlarging the dual lattice (of excitations) and of introducing a selection rule (the 
Znjj pairing rule). In fact, there exists a decomposable sublattice L c F of the same 
dimension spanned by the vectors 

{e\a.} , e^ = nniq - oj) = Qdn, (17) 



where are the basis vectors of the neutral sublattice F^y c F, i.e., the simple 
roots of the Cartan lattice. This sublattice indeed splits into 2 mutually orthogonal 
sublattices: 

L = ZQdH®Tw . (18) 
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We have the inclusions L c F c F* c L* and the physical lattice F can be decom- 
posed in terms of its sublattice L as follows 



F = |7 = A + 5£^j AgL, 0<5<n//-l|, 



(19) 



i.e., Vy e F, 3! i e L, and 3! ^ e such that y = A + sq^^. Note that L* /T* ^ 
T/L ^ Z„jj since the determinants of the Gram matrices of L and F (which give the 
number of sectors of the corresponding rational conformal theory) are: 

1^1 = dl\Qf\Cw\' = nndn det(C^y) = 4|F|. (20) 

The spin-charge separation of excitations is achieved in the decomposable dual lat- 
tice L*, whose physical points (corresponding to the points of F*) obey a " mod hh"- 
selection rule: if A G L*/L then A is a physical excitation (i.e., A G F*/F) iff 

\fAeL*/L: Aer/rcL*/L « (^^J^) G Z. 

In other words, if we write the excitation charge A in the basis of L* 

^ = /e;+A, A = 2a'A,., I, a' £ Z (21) 



with A. being the fundamental weights of the Witt sublattice, then A is a physical 
excitation if 

(^i'-) ^ ~~ ('^"{oA^ - /) G Z <^ nnioAA) = I mod Hh, (22) 



where co = 'Zji ^'^i with w' e Zis the neutral component of the electron's topolog- 
ical charge (14). This condition reduces the number of elements in L* /Lhy a factor 
of njj to that of F*/F which is the number of physically different quasiparticles. 
To summarize, we first extend the dual lattice F* to L* obtaining the spin-charge 
separation in the latter and then identify the physical points in the extended lattice 
L* by a Z„^ pairing rule. 



N 

2. 6 The CFT characters for the lattice extensions ofu(l) 



Because of Eq. (19) the Hilbert spaces H^of the original F*/F model can be written 
as direct sums of the Hilbert spaces 'H^ of its decomposable sublattice (18) 

- yer - - s mod hh ^ei - j mod hh - - 2ei 
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N 

where 'Vy* is the m(1) Verma module (A^ G N is the rank of the lattice F) over the 
lowest-weight state I7*). Therefore the F-characters would have a similar decom- 
position into L-characters 

7*+%] s mod hh 

where the factor n^f in front of ^ appeared because 2 = <?* in F* but in the decom- 
posable dual lattice it is 2 = nz/Cj e L*. Because the sublattice L is decomposable 
the characters of the "H^ modules are simply products of those for the correspond- 
ing charged and neutral sectors. Writing explicitly y* = /e^ + A we have 

Xi,k{t,C^T)= 2 Ki^,dHi'^,nH^;nHdH)da^*h{T) (23) 

s mod Hh 

where Ki{t, Hh^; nndn) is the character of the (charged) m(1) sector, i.e., the rational 
torus partition function [15] 

^ 00 

Ki{T, ^■m)=—-y ^5("^^)'e^'^'-^("^^), //(r) = f[ (1 - ^") , (24) 

^^^^ n€Z n=l 

(77 is the Dedekind function [15]) while chA(T) is that of the neutral sector of the 
decomposable lattice L and A satisfies the Z„^ pairing rule nff(ct>|A) = / mod hh- 
Note that the electron charge vector (14) acts on the neutral sector by the s-i\\ power 
of the simple current io which in an abelian CFT is taking the form o/*A_ = sco + K. 
The characters expressions (23) are valid for the FQH state corresponding to the 
arbitrary topological charge lattice F. 

2. 7 Neutral reductions and stability criteria 

The classification of the FQH universality classes in terms of chiral quantum Hall 
lattices is motivated by the rationality of the filling factors and of the electric charge 
spectrum. In addition, according to Eq. (16) the neutral component io of the elec- 
tron's topological charge should also have a rational CFT dimension depending 
on the filling factor (see Eq. (13)). That is why the electric properties of the FQH 

N 

liquid are well-described by the topological charge lattice extension of the m(1) 
CFT introduced in Sect. 2.3. However the electric properties, which are basically 
what is measured in the experiment, are not enough to fix a unique rational CFT 
corresponding to a given FQH universality class. For example, we can start with 

the lattice extension of the m(1) CFT and^oject some degrees of freedom in the 
neutral sector keeping the same electric m(1) part. This would certainly preserve 
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the electric properties and produce another candidate CFT for the same univer- 
sality class. Two important consistency conditions are that the neutral projection 
preserves the neutral component oj of the electron topological charge vector, for 
locality reasons, and the pairing rule which is a general property of the FQH 
excitations due to the spin-charge entanglement. The advantage of the neutral pro- 
jections is that they generically reduce the total central charge and therefore may 
produce better candidates for the same universality class according to the stability 
criteria formulated in Ref. [11]. The neutral projections also preserve the structure 
of the characters (23) only the neutral characters ch^^^A have to be replaced by the 
corresponding projected ones. This is important because the neutral characters give 
non-trivial contribution to all thermodynamic quantities. An example of a neutral 
projection which we shall describe in Sect. 3 is the affine coset reductions which 
realizes the parafermion FQH states in the second Landau level [8]. 

The FQH effect phenomenology seems to be in favor of the neutral projections: 
the main quantum numbers expected in single-layer disk QH samples of polarized 
electrons are electric charge and orbital momentum, while in the lattice extensions 



of m(1) there are A'^ separately conserved charges. The projection would remove 
the superfluous charges keeping only the electric charge and the orbital momentum, 
and in some cases would preserve discrete quantum numbers related to W„ algebras, 
like in the coset projection in Sect. 3, due to the Z„^ pairing rule. 



2.8 The Aharonov-Bohmflux: twisting the electron 

Let us introduce the AB magnetic field piercing the FQH disk sample at the origin 
in the z direction which is produced by a vector potential of the form 



is the length of the radius vector r = {x,y) in the plane, is the unit vector along 
the direction of the polar angle (p and <^ e R is the AB magnetic flux in units of the 
flux quantum <^o = hje. In this section we shall see how the partition function is 
modified in presence of AB flux. 

First, because rot A = ^ 6^^^{y) is zero everywhere except for the origin, the electron 
field iff^iiz) in presence of AB flux can be written as 



N 



A = -- — - {-y, x) = -- — e^, where r = ^f^Ty^ 




(25) 



z 



il4i(z) = exp -j- 




.dr hAei(z) = z"Vei(z), 0eM, 



(26) 



* 
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(e = h = 1) where i//s\(z) is the electron field operator in the absence of AB magnetic 
field, provided that the path of the line integral does not encircle the origin. In 
the computation of the line integral in Eq. (26) we have used z = rexp(i(p) and 
have chosen a path along the cp coordinate line ^ , i.e., dr = re^dcp, ending at 
(p = arg(z), the polar angle of the coordinate z. The second equality in Eq. (26) is 
a simple special case of a more general procedure called twisting [18,19]. The AB 
flux modifies only the m(1) part of the electron operator (15) which, for any QH 
system (integer and fractional), is realized as the m(1) vertex exponent defined by 

^<^t^>(,} fJ. e'"^^"fe), /0(^')(^) = ± y jj—, (27) 

^ n 



where Ua are outer automorphisms of the m(1) current algebra satisfying Ua Up = 
Ua+p for a,p G R, C/o = 1, {Uj = and 

[i„, Ua] = aUa 5„fl. (28) 



Therefore we can use the well known results [18] for the twisting of the normalized 
m(1) current, u{l) Sugawara stress tensor and normal ordered exponentials. The 
twisting with a twist yS is the transformation [18] Tr^g : A ^ UpA U-p for any 
operator A from the chiral (observable) algebra; in particular twisting of the vertex 
exponents gives 

.g,-a^<^'a). j!^ Up:e"'^'''^''^-.u-p = z-"^-.e"''^'''^''^-. (29) 



which can be computed directly by using the definition (27) and the commutation 
relations (28) and (7). Because the twisted electron operator 

np ((Aei) (z) = Tip (:e-^^'''^^^) 0^(z) = z-^^ijf.iiz) (30) 



has to reproduce Eq. (26) the value of the twist p corresponding to the addition of 
AB flux can be determined by 




Then the normalized m(1) current (7) is modified as follows [18] 

J{z) ^ UpJ(z)U_p = J(z)-^, (31) 

z 

^ in the CFT we conventionally choose r = 1, however, the results can be extended to 
arbitrary r by analytic continuation of the CFT correlation functions 
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which again can be computed from the definition (7) of the current and the com- 
mutation relations (28). Next, the twist of the electric m(1) current (8) is 

r\z) ^ r\z) + ^ « Jf ^ Jt + VHcp6n,o. (32) 



Furthermore the twist of the stress tensor T^'^'iz) of the m(1) part, defined by the 
Sugawara formula (9), can be derived from Eq. (31) to be [18] 

z 2z^ 



so that the total stress tensor is modified as follows 

2 2 

T(z) ^ T(z) + ^r\z) + ^ ^ L„ ^ L, + cpjf + y^^c5„,o.(33) 
z 2z^ 2 

Remark 1 The twisted electron 7r^(iAei) has a different CFT dimension with respect 
to the untwisted CFT generator Lq, i.e., different spin and statistics so that it is 
not local with respect to Lq and therefore their correlation functions would not 
be single-valued in the electron coordinates. However, with respect to the twisted 
generator n^(Lo) = A) + ^Jq + Vh'P^I'2 the twisted electron field n^(i/fei) is local and 
has the same CFT dimension like the untwisted one with respect to Lq. 

Using Eqs. (33) and (32) for the zero modes Jf^ and Lq of the twisted electric current 
and the twisted stress tensor respectively and putting them into the partition func- 
tion (4) we get the following expression for the CFT partition function in presence 
of ABflux<^ 



Z^(t, = tr exp !^2ni t |lo + (PJ^ + vh^ " ^) + ^ (-^o + ^//<^) 
. (inivH i^T + J) Z(t, ^ + (f>T). 



= exp 



(34) 



The simple expression (34) for the partition function with AB flux is valid for any 
quantum Hall state and can be directly used for the computation of the free energy 

(1) in presence of flux, hence, of the persistent currents, AB magnetizations and 
magnetic susceptibilities for temperatures small on the scale of the corresponding 
energy gaps. Thus Eq. (34) which is a central result of this paper is the cornerstone 
of the eff'ective CFT description of general mesoscopic phenomena. 
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2. 9 Laughlin 's spectral flow and Cappelli-Zemba factors 



Two of the assumptions in the CFT approach to the FQH effect, based on the FQH 

phenomenology, are that there is a minimal electric charge in the quasiparticle spec- 
trum depending on the filling factor and there are only a finite number of topologi- 
cally distinct quasiparticles [2]. This means that the CFT describing the FQH states 
should be rational. While this has not been rigorously proven it seems reasonable 
and most of the FQH states for which the CFT are are well established are indeed 
rational. The only exception are the non-rational W\+oo minimal models [20] for 
the Jain series of FQH states. For the CFT description of mesoscopic phenomena 
in FQH systems rationality is not necessary, however another condition, called the 
V invariance representing the Laughlin spectral flow, is crucial in this case and we 
shall describe it in this subsection. 

The rational CFT characters (5) belong to a finite dimensional representation of 
the "fermionic" projective modular group [8,11,14] r(2) c PSL(2,Z), i.e., they 
should be S and covariant [14] where S and are represented by the following 
transformations of the modular parameters 



In addition, due to the specifics of the FQH effect there are two more covariance 
conditions to be satisfied by the characters [14] which represent certain electric 
properties of the FQH liquid: the annulus partition function must be also invariant 
under the U transformation of the modular parameters [14] 



The V transformation (35) can be interpreted as adding adiabatically one flux quan- 
tum to the FQH system as can be seen from Eq. (34) if we put 0=1. Therefore 
Eq. (35) represents the Laughlin spectral flow [21,14], i.e., the mapping between 
orbitals, corresponding to increasing their orbital momentum by 1, when threading 
adiabatically the sample with one quantum of AB flux. This mapping preserves the 
entire quantum spectrum of FQH excitations and should be an exact symmetry of 
the partition function [21,14]. In a FQH state, as a result of this flow a fractional 
amount of electron charge is transferred between the two edges of the annulus, and 
this is the Laughlin argument describing the Hall current [21]. However, as noted 
by Cappelli and Zemba, the CFT partition function for the annulus FQHE sample 




U: ^^^+1 



and under the V transformation 



(35) 
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is not V-invariant alone [14] because the transformation (35) changes the absolute 
values of the CFT characters Xa in Eq. (4). To preserve the norm of the charac- 
ters corresponding to a general FQH state, and restore their V covariance, these 
authors introduce special universal non-holomorphic exponential factors 

exp I -TTVh I , (36) 



multiplying the characters (5), which we call the Cappelli-Zemba (CZ) factors. 
The chiral partition function (4) also becomes V invariant only after multiplying 
the characters with the CZ factors (provided Re r = Re ^ = as appropriate for the 
chiral partition function) 

Z(T, = e-'^''-^ Z(T, =^ Z(T, ^ + r) = Z(T, 0- (37) 

In order to prove the V-invariance of Z(t, ^) we first use the following property of 
the ^-functions 

e-^«^^Ki(T, nni^ + (frr); nHdn)) = Q-'^"^ Ki,„,^(t, HhC, nHdn), (38) 



which follows directly from the definition (24) of the K functions and CZ factors 
(36), applied for ^ = 1. Second, because of the pairing rule (22) it can be shown 
that the labels of the neutral characters A = are among the fundamental 

weights of the Witt sublattice so that ch^^^A = ch/ mod nH,s and therefore chi+n„,s = 
chi^s which proves the "quasiparticle's spectral flow" 



Finally the V-invariance of partition function (4) follows from the invariance of the 
sum under the shift of the m(1) label / (corresponding to a shift of A in Eq. (4) ). 

The CZ factors were originally interpreted [14] as adding constant capacitive en- 
ergy to both edges which makes the ground state energy independent of the edge 
potential Vq. For the single-edge disk sample, however, it is more intuitive to inter- 
pret the CZ factors directly in terms of the AB flux encircled by the edge. In this 
paper we are going to show that these factors play much more important role in 
the equilibrium thermodynamic phenomena, such as the persistent currents and AB 
magnetization. 

Now we can show that the partition function with AB flux (34) is invariant under the 
Laughlin spectral flow and is essentially equivalent to the procedure of Ref. [14] of 
multiplying the conventional partition function by the CZ factor. First, because the 
partition functions (34) and (4) are chiral (corresponding to a single-edge disk FQH 
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sample) the modular parameters r and ^ are purely imaginary, i.e.. Re r = Re ^ = 0. 
Then the exponential factor in Eq. (34) can be written equivalently as 



and we note that exp ^-;rvij Im t + ^) j = exp|-;rv^f Imr^^^) j for ^' = 
^ + (pTSO that the partition function with AB flux (34) is simply 



Z^(T, = Q"""" "^^^> Z(t, ^ + 0T) (40) 



where the exponential prefactor is now (p-independent and Z(t,^) is defined in 
Eq. (37). Next, for integer (p, the invariance of Z^(t, ^) under (p —> ^ + I follows 
from Eq. (37): 



Z^+i(t, = e'^«^^(^)' Z(t, ^ + <pT + T) ^=^^ e"^»M^f Z(t, ^ + (pr) = 

This result is an illustration of the Bloch-Byers-Yang theorem [22] which states 
that the partition function of any non-simply connected conducting system threaded 
by magnetic flux is a periodic function of the flux with period of 1 flux quantum. As 
can be seen from Eq. (39) the spectral flow (35) maps one character into another and 
in order for the partition function to be V-invariant it must include all inequivalent 
sectors which are connected by the action of the spectral flow Eq. (35). 

To summarize, the arbitrary flux (p threading procedure within the rational CFT 
framework, is represented by the following transformation of the parameters in the 
partition function 

^^^ + 0T, 0eR (41) 



provided that the CZ factors have been included. Then the partition function for 
the edge states of a FQH system in presence of AB flux represented by the vector 
potential (25) is given by either Eq. (40) or Eq. (34). 

Finally we can make the following identification of the modular parameters in the 
context of a disk FQH threaded by AB flux: as we have already seen Re r = Re = 
0, and Imr ~ llksT where the proportionality constant sets the temperature scale. 
We choose the scale Tq from Ref. [4] so that 

Tn hVp 

T = in^, ^ = cPt, 7-0 = —^, (42) 
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where T is the absolute temperature, L = 2nR is the circumference of the edge, 
kg the Boltzmann constant and (p the magnetic flux (in units of h/e = 1) threading 
the FQH disk. Note that for ^ = + 0t, which will be our case, the two partition 
functions, Eq. (34) and Eq. (37) coincide, i.e., Z^(t, 0) = Z(r, (pr) so that both can 
be used for the computation of the persistent currents. Let us note that Eq. (41) 
together with the CZ factors (36) have already been used in Refs. [6] and [7] for 
the computation of the persistent currents, AB magnetization and susceptibility. 

2.10 Low-temperature asymptotics of the partition function 

The low-temperature limit T ^ in the partition function corresponds to the trivial 
limit ^ ^ 0. Therefore we can keep only the first terms in q in the linear partition 
function (4) which come from the vacuum character and from those with the min- 
imal non-trivial CFT dimension. In this paper we shall assume that the minimal 
nonzero CFT dimensions come from the sectors containing either one quasiparticle 
(qp) or one quasihole (qh) but not both. In this case the partition function can be 
approximated by 

^(J,0 ^=^^ Q-'''"^'^ (Ko(T,nHC,nHdH) choir) + 

+ Ki(t, nn^; nHdn) chqh(T) + K-i(t, nu^; nudu) chqp(T)), (43) 

where hh and dn are the numerator and denominator of the filling factor (10), re- 
spectively. This assumption is natural for the FQH states for which the CFT dimen- 
sion of the quasihole/quasiparticle, which carry the minimal absolute value of the 
electric charge (corresponding to / = +1 in the A" functions) in the spectrum, is the 
minimal one allowed in the CFT. Notice that the maximally symmetric quantum 
Hall lattice CFT for the Jain series labelled by {2p + l|'A„_i) does not satisfy this 
assumption, however this is not a problem since it is almost certain that this CFT is 
not the correct one [5]. In Eq. (43) cho(T) is the CFT character of the vacuum rep- 
resentation of the neutral sector while chqh(T) and chqp(T) are those of the neutral 
sectors containing one quasihole and one quasiparticle, respectively and we note 
that only the leading powers of q in the characters are important for T — > 0. 

To include the AB flux we apply the arbitrary flux transformation (41) to the par- 
tition function (43) and use the property (38) of the ^-functions (24) for ^ = 0. 
Moving the flux dependence into the indices of the J^-functions is technically more 
convenient since it guarantees the reality of the partition functions when computed 
numerically for arbitrary (p. Keeping only the first terms in Eq. (43) and skipping 
the 0-independent terms ^ containing the 77-function and the central charge we can 

^ for the computation of the the persistent currents we intend to differentiate the logarithm 
of the partition function with respect to according to Eqs. (1) and (2) so that the 0- 
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write ^ 



Z(T, (P) ^ =^ q-^'^' (l + (^Gqh^ + ^-Qqi.^)} = 

= e-'^^^^^'jl +2e"S^cosh|2;r2y!2qh0)| (44) 

where 

iqh = 2;r2^BroAqh, Aqh = + Aj^^J, (2qh = :^. (45) 



In the above equation Sqh is the proper quasihole energy, Aqh is the total CFT dimen- 
sion of the quasihole (with A^^^^ being its neutral component) and Q^^ is the quasi- 
holes electric charge. In the derivation of Eq. (44) we have also used cho(T) - q^, 

AO) 

chqp(T) = chqh(T) ^ ^ ii- for ^ — > (again dropping the central charge term). 

Remark 2 The partition function (44) is finite in this approximation since for all <p 
one has v//0^/2 + Qqh(f> + Aqh > because the discriminant of the above quadratic 
equation is -Ivh^^^I — ^ provided that h^^l > 0. However, when in the low- 
temperature limit we take the logarithm in Eq. (1) we are going to use the ap- 
proximate expansion ln(l + Z) ^ X which is valid for X I. In our case, for 
r «: To 

Z = 2e-^ cosh(-2;r2^eqh0) « 1 iff |0| < ^ and Aqh - ^Qqh > 0. 



The last condition is always fulfilled in the FQH states, such as the Laughlin states 
and the Z^t parafermion states which satisfy the quasihole charge-statistics relation 
[13] 

2Aqh = !2qh. (46) 



As we shall see in Sect. 4.3 the low-temperature asymptotics of the persistent cur- 
rent's amplitude derived from Eq. (44) is logarithmic 

Im^T) = -Vh—<1-— 1+ln — — 
r«ro 2 L [ gqh [ ynnkBTj 

independent terms would not contribute 

^ after including the flux the indices of the K functions are modified so that e.g., K2+nH(p 
may dominate over K-\+nH<P for -1/2 < < in the T ^ limit. However, in this case 
both functions can be neglected when compared to Ki+„jj^ since the latter gives the smallest 
minimal CFT dimension 
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and shows that the persistent current decays logarithmically with increasing tem- 
perature due to the probability for thermal activation of quasiparticle-quasihole 
pairs ^ which reduces the radial electric field that is responsible for the appearance 
of the azimuthal persistent current. 



2.11 High-temperature asymptotics: modular S matrix and quantum dimensions 
of quasiparticles 



The high-temperature limit 7 ^ oo is rather non-trivial since the modular param- 
eter q = exp(-2n^To/T) 1 is at the border of the convergence interval for the 
partition functions and therefore cannot be taken directly. However, one could use 
the advantage that the complete chiral partition function is constructed as a sum of 
RCFT characters, so that 

N 



where the modular parameters (r, ^) and (r', are related by 



S : 



r = 



V 

t' 



t' = -'-= i^ 



(48) 



Apparently Eq. (48) connects the high-temperature and low-temperature limits of 
the partition function. Note that the usual common phase factor [14] in front of the 
S matrix is trivial since exp (/Trv^fRe (^'^/r')) = 1. Now the new modular parameter 
vanishes in the high-temperature limit 

^' = e^^^'^' = exp|-^j ^ when T ^ oo (49) 



so that it would be enough to keep only the leading terms. After the S transforma- 
tion (48) the partition function (4) becomes ^ 

N N 

Z{T, = Yj P^Xa{t', n Pa = J]Su'. (50) 

^=1 A'=l 



note that the proper quasihole energy appears with a factor of 2 
^ the CZ factors (36) are trivial after S transformation exp(-;rv^f(Im^')^/IniT') = 1 since 
C ^<f>, i.e.,Im^' = 
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Obviously the universal quantities Fx are supposed to play an important role in the 
high-temperature analysis. 

Remark 3 Note that Fx are always real. The reason is that the S matrix is sym- 
metric Sx,p = S and satisfies Sx,p- = {s p,xj [15] where p* is the the label of 
the charge conjugated representation of\p), i.e., |p*) = C|p) with C = S^. There- 
fore, for self-conjugated weights, p* = p, the S matrix has real rows and columns 
Sa,p = iSx,p) (for all A) while for p* = a p the rows and columns corresponding 

to p and (T are complex conjugated S xxr = a,p) (for '^11 Because the rep- 
resentations are either self-conjugate or enter in conjugated pairs the sum of the 
elements of an entire row or column, i.e., over over all representations, will always 
be real. 

Since the complete CFT contains a m(1) factor it follows from the properties of the 
simple currents that 

Fxi^O « QM = 0. 



In particular Fi, where 1 denotes the vacuum representation, is the following sum 
Fi=SnJ]Dx, where Dx = > 



are the quantum dimensions of the topological excitations of the FQH fluid. Note 
that the quantum dimensions have several nice properties which justify the name 
"dimensions" 



where A^'^v'' are the fusion coefficients [15]. The second of the above equalities, 
which follows from the unitarity of the S matrix, shows that the vacuum in a ra- 
tional CFT does know about all topologically distinct representations because the 
matrix element S n depends only on the vacuum module. 

Using Eq. (50) we are going to show in Sect. 4.4 that the high-temperature asymp- 
totics of the persistent current's amplitude in general has the form 
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with a universal exponent a which can be written in the form 
Vh 

where is the neutral weight with the minimal CFT dimension satisfying the 
pairing rule P[A.J = dn mod hh, see Eq. (78). According to Eq. (51) at high tem- 
peratures the mechanism suppressing the persistent currents is different from the 
low-temperature one: the amplitudes of these currents decay exponentially due to 
the effects of thermal smearing and decoherence. 



3 Application: the parafermion FQH states 



3.1 The CFT partition function for the It^ parafermion FQH states 



Following Refs. [8,13] we denote the CFT for the parafermion FQH states (PF^t) 
in the second Landau level by [8,13] 

yH = ^ (52) 



where vh is the filling factor in the last occupied Landau level. This CFT contains 
the standard m(1) factor describing the charge/flux degrees of freedom and a neutral 
Zyt parafermion sector which is realized as a diagonal affuie coset [8] 



su{k)i e su(k)i (0) lik - 1) 

= — ' ''p^ = TT2" ^^^^ 

where Cpp is the central charge of the coset. The affine coset realization of the 
parafermions has a long history, please see Refs. [23-25] and the references in [8]. 
The primary fields of the coset (53) are labelled by [8] 

Ap+Ap> 0<fi<p<k-\, (54) 



where A^, for /i = 0, . . . , ^ - 1 are the suik) fundamental weights (Aq = 0). Because 
the numerator of the filling factor is k the two factors in Eq. (52) are subjected to 
a 2,k pairing rule [8], in agreement with the general discussion in Sect. 2.5, i.e., an 
excitation 

0(/; A^ + A^) = :e'^*^"(^^ ^^^{h.^ + A)(z), 0<fi<p<k-l (55) 
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of the PFjt model with labels (/, + A^), where / is the m( 1 ) label and O^^ (a^ + A^) 
is a Z^-parafermion (primary) field [8,15], is legal only if 

P [O^^ (a^ + A,)] = / mod k, where P [^^^ (a^ + A,)] = + P mod k 



is the parafermion 'Lk charge. The CFT dimensions of the coset primary fields 
(a^ + Ap) 

have been found in Ref. [8] in a compact form 



The independent characters xi,p of the chiral CFT for the parafermion FQH states 
are labelled [8] by a pair of integer numbers (/,p) restricted by 

/mod/c + 2 and p = Q,...,k-\, where p>/-pmod^ (57) 



which determines the topological order to be {k + l){k + 2)/2. The chiral partition 
function with the CZ factors (36) sitting explicitly in front of it is 

Zk{T,0 = ^-'^'''^ Yj Z (58) 

/ mod k+2 p>l-p mod k 



where 

k-l 

XiAr, = J] Ki^sik^T, k^; k(k + 2)) ch (a^_^^^ + A,^ J (r) (59) 

i=0 



The characters Ki(T,k^;k(k + 2)) of the m(1) sector are defined by Eq. (24). The 
neutral characters ch(A^ + A^) which are the characters of the coset (53) can be 
found using the general approach [24] are labeled by the weights (54) and fol- 
lowing Ref. [8] can be expressed in terms of the so called universal chiral partition 
functions (see Ref. [8] and the references therein) ch^ after the identification Q = p, 

°° /v«-C"'-(n-A<^) 

che(T) = /-^ V -\ Q<a<Q. (60) 

[n\=Q mod k 



In the above equation, n = (ni, ^2, . . . , rik-i) is a ^ - 1 component vector with non- 
negative entries, [n] = 2m) * is its A^-ality, A^- is the CFT dimension (56) of the 
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PF^: primary field labeled by[8] Aq + A^, cpp = Cpp is the parafermion central charge 
from Eq. (53), 

n 
7=1 

and is the inverse of the su(k) Cartan matrix. The partition functions with AB 
flux, is computed according to Sect. 2.8 by applying the flux transformation (41) 
to the partition function (58)and using again the property (38) of the i^-functions 
for ^ = 0. In Sect. 4.2 we shall use the CFT partition function (58) with AB flux 
to compute numerically the persistent currents in the parafermion FQH states for a 
wide range of temperatures as well as to find analytical formulas for the low- and 
high- temperature limits. 



3.2 The modular S matrix for the Zi^-parafermion FQH states 

For the analysis of the high-temperature asymptotics of the persistent currents in 
the parafermion FQH states we shall explicitly need the modular S matrix for the 
corresponding CFT. The computation of the S matrix for the m(1) x PF^- model is 
a bit complicated by the fact that the characters (59) are not simply products of 
characters of its neutral and charged sectors but sums of these. Nevertheless the S 
matrix can be found by carefully using the properties of the simple currents [15,17] 
and following the standard Goddard-Kent-Olive approach [26,15] for the coset. In 
order to compute the full S matrix we note that the full PF^ characters (59) can be 
written in terms of the simple current 

J = :exp (/-A±|^0(')(z)]: ® ^^^(A^ + A^)(z), (61) 



(coinciding with the physical hole operator) which acts on the primary fields (55) 
according to 

J * {-..^imT/''^^^ ^ 0-(A^ + Ap(z)) = -..m/'^^^ ^ 0-(A^,, + A^^.Kz), 

as the sum over the orbit of the product Ki ch(A,_^ + A^), i.e., 

k-i 

Xi,piT, = YjJ'* {Ki (t, k^; kik + 2)) ch(A,_^ + A,)(t)) . 

4=0 
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Next we use the general property of the simple currents' action over the S matrices 
[15,17]: 

Sj.AM = e-^'^'G^^^^^A.M, Qj(M) = A{J * M) - A(M) - A(J^) mod Z (62) 



where A and M are the labels of any two coset primary fields, A(A) is the total CFT 
dimension of the field labeled by A and Qj(M) is the monodromy charge [17]. 
Because ^ is a simple current the fusion product J'* M contains only one primary 
field so that the monodromy charge is well defined. Then we apply the modular S 
transformation to a single product of characters Ki ch(A^_^ + A^), namely the first 
term 5 = in the character (59) which gives 



SKi(T, H; kik + 2)) ch(A^_^ + A^Xr) = ^''''i'^ ) J] S\]I,K,{t, k^; k(k + 2)) 

/' mod k(k+2) 

X 2 ^I;rch(A^ + A)(T), 

0<ti<v<k-l 



where 



S = ^ exp I -2;rj — | , with /, /' mod k{k + 2), (63) 

ik(k + 2) \ k(k + 2)) 



is the S matrix for the u(l)j^(^j^+2) ^ n,v the one for the diagonal coset (53) which 
will be discussed in more detail in Appendix B. Next we apply the S transformation 
to the full character (59) using the property (62) of the simple currents (61) to get 



SXi, = ^"'"^'^^^f, Z 2 exp[-2;r/(GlJkr) + Q^.(A +A))]x 

s=0 I' mod k(k+2) 0</i<v<fc-l 

X 5;;^ s'^^'^KKt, kC, k(k + 2)) ch (a^ + aJ , (64) 
where the monodromy charge Q'^^^ for the m(1) part is (using A'^^^(/) = P/2k(k + 2) ) 

!2y.(/') = A(i> (/' + s(k + 2)) - A(i) (/') - A('> (s(k + 2)) = ^, 



and that for the coset part is (using Eq. (56) for the CFT dimension in the coset) 
S.. = A" (V. - A J - A" (A, . A J - A" (A, . A.) = . 
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Because the only 5-dependence is under the exponent of the monodromy charges 
the sum over s can be taken expUcitly and gives a J-function 



^ exp |-27r/ j = k6(l' = fi + v mod k) . 



Plugging this (5-function back in Eq. (64) and splitting the summation variable 
/' mod k(k + 2) into 

r = /" + s"(k + 2), where s" mod k and /" mod k + 2, 
we get (setting v =: p' ^ ix = I' - p' mod k < p') 

^Xl,p-^ / J / J ^ U"+s"(k+2) ^ /"-p'+/'(^+2),p' ^ 

/" mod *:+2 i=0 p'=0 

X is:/v+,,,(4,+2)(T, k^, k(k + 2)) ch (A,„+,„(i+2)-p' + Ap') • 

In the above equation we can change s" -p' =: p" and interpret the sum over s" as 
a simple current's orbit because 

A,„^ „„^o^ , + A , = A,„ „^ „ + A „^ „ = * (A,„ „ + A „ ) . 

— l"+s"(K+2)-p' — p' — l"-p"+s" — p"+s" \ — l"-p" — p" I 

Using once again the property (62) of the 5 matrix under the action of the simple 
currents, i.e., 

^ l/"+."(^+2) = exp -27ri— S J/, and S i„_p„^,,^p„ = exp -^ni—— S ,„_^„^^ 



one could see that the two exponential factors cancel each other and the sum over 
s" exactly reproduces the character ;^:'///_p// according to Eq. (59). Thus we finally get 
the following compact expression for the total S matrix of the m(1) x PF<: model 

k-\ ^ 
Sxi,p= J] ^Sr,p"Xi",p" with s\^p„=kS^^„Siffp„^p,„ (65) 

/" mod k+2 p"=0 



where I - p mod k < p and the prime over the first sum means that the summation 
over /" and p" is restricted by (/" -p" mod k) < p"; the m(1) modular matrix S '-^^ is 
given in Eq. (63) while the coset S matrix can be written explicitly in the basis of 
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coset primary fields ^ (A^ + A^) labelled by + A^ with 0<//<v<^-las 



0<H<v<k-\ 
Q<p<a <k-\ 



(66) 



In Eq. (66) S^-^^^, is the S matrix for su{k)2 which can be computed as a sum over 
the su{k) Weyl group [27] 



:k(k-l)/2 ^ 



exp 



|a + p w(A' + p)) 



-27n" 



A: + 2 



(67) 



The derivation of diagonal coset S matrix (66) following the Goddard-Kent-Olive 
construction [26,15] is sketched in Appendix B. Eq. (66) shows that the S matrix 
for the diagonal coset (53) is related to that of the su(k)2 current algebra, by charge 
conjugation and the action of simple currents which are automorphisms of the fu- 
sion rules, i.e., the two S matrices are equivalent. Therefore not only the irreducible 
representations of the coset (53) and su(k)2 are labelled by the same weights (54) 
but the fusion rules of both CFTs are identical. This explains why the quantum 
group structure of the affine coset (53) is the same as that for su{k)2 (see Ref. [28] 
and the references therein), i.e., 

I n \ (t - q~"^ 

Uq {sl{k)) with q = exp l^lg = i _ _i • 



4 Persistent currents 



4.1 An intuitive picture 



Consider the Corbino ring sample shown in Fig. 1 threaded by AB flux that is 
changed adiabatically at zero temperature ^ . We shall assume that the huge back- 
ground homogeneous magnetic field Bq = -Bq e^ (perpendicular to the plane), which 
points in the negative z direction, has driven the electron system into an incompress- 
ible FQH state on the ring so that cTxx = and cr^y = vne^/h. The actual magnetic 
field is measured with respect to 5o, i.e., B = Btot - Bq. When no flux is threading 
the sample, i.e., 5 = 0, the currents are concentrated on the edges and have the 
same magnitude but opposite directions so that the net current is zero. For a sam- 
ple with negative-charge carriers the (non-mesoscopic) current on the inner edge is 

^ this process is reversible on each step 
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Fig. 1. The two-dimensional Corbino disk in ring geometry (CD-R) and the persistent cur- 
rents on both edges 

counter-clockwise while that on the outer edge is clockwise as shown in Fig. (1). 
An intuitive way to determine these directions is provided by the skipping orbits 
semi-classical picture: the electrons in the plane move along clockwise circles be- 
cause the magnetic field points downwards. In the bulk the electron orbits are fixed 
complete circles so that there is no electron drift, hence no net current. However 
on the inner edge the circular orbits are cut and elastically reflected from the edge 
becoming open curves resulting in a clockwise drift of the inner-edge electrons pro- 
ducing a counter-clockwise current /_. On the outer edge the drift of the electrons 
in the cut orbits is counter-clockwise so that the outer edge current /+ is clockwise. 
Quantum-mechanically the situation is quite similar only the eff"ect of the edges is 
represented by a confining potential. Note that the role of the confining potential is 
fairly important - without it there would be no mesoscopic eff"ects. 



Introducing adiabatically the AB flux creates by Faraday induction an azimuthal 
electric field E = -t--^ e^.. which cannot create an azimuthal bulk current since 

Inr at V 

a XX = but induces a radial current pulse due to the non-zero cr^y 



1 = Vh 



1 
-1 



# 1 



Basically the radial current pulse creates a radial polarization which on its own in- 
duces an axial quasi-permanent electric current which is called persistent since it 
is observed for more than 10^ s in real samples. In more detail, the current pulse i. 
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corresponding to the Hall current which is purely bulk in this geometry, transfers 
electric charge = -evH^4> between the two edges which destroys the perfect 
balance between them. This pulse increases the carrier density on the outer edge 
while decreasing that on the inner one and therefore increases the outer-edge cur- 
rent while decreasing the inner-edge one by the same amount. As a result there 
appears a net current which is diamagnetic, i.e., clockwise, for small positive A0, 
that is in the direction of the electric field (under the reasonable assumption that 
the Fermi velocities do not change with the AB flux). Note that the mesoscopic 
net current fluctuates with period 1 and can therefore be also opposite to the elec- 
tric field as the flux changes. The above argument also implies that the mesoscopic 
persistent currents (i.e., the changes in the edge currents due to the additional AB 
flux) on both edges have the same direction (clockwise) and both produce negative 
magnetic moments. 

The typical form of the mesoscopic edge currents for 7 = can be obtained from 
the charge transferred between the edges: 

A2 evp 
AI = eVfAn = Vf — — = ±vh— — (p. 



where we have divided AQ by the circumference L in order to get the one-dimensional 
charge density on the edge (e L An = AQ). Thus we get 



(out) (in) 

evp evp 



We note that the mesoscopic edge currents are just small fluctuation to be added 
on top of the much bigger non-mesoscopic edge currents [9] produced by the huge 
background magnetic field. 



4.2 The thermodynamic CFT derivation 



According to Eq. (2) and the discussion in Sect. 2.8, the equilibrium chiral persis- 
tent current in the FQH system can be computed directly from the CFT partition 
function by 

I = i-\kBT-^\r\Z{T,(bT), where T = in—, Tn = , (68) 

Vf is the Fermi velocity and L = InR the circumference of the edge. Because of 
Theorem 3 in [22] the partition function of the system should be an even periodic 
function of the flux with period 1. Therefore the persistent current (68) must be an 
odd function with period 1 so that it would be enough to consider the half interval 
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-1/2 < < and continue it by 7(0) = -I(-(p). The V-invariance of the CFT parti- 
tion function (58) for the PFyt states implies that the corresponding chiral persistent 
current should be periodic in <p with period at most 1. Any other value, such as 1 /2, 
1/3 etc., of the persistent current's period which is smaller than 1 is a signal for 
some broken continuous symmetry which is not allowed in unitary field theories in 
1 + 1 dimensions [29,30]. Therefore we expect that the period should be exactly 
one flux quantum. The plots of the chiral persistent currents in the first several PF 
states computed numerically from Eq. (68) for -1/2 < (p < 1/2 at temperature 
T/Tq = 0.1, given on Fig. 2, indicate that all currents are periodic in with period 

Fig. 2. Persistent currents in the A: - 2, 3 and 4 parafermion FQH states, as functions of 
the magnetic flux within one period, computed numerically for T/Tq = 0.1. The flux is 
measured in units h/e and the current's unit is evp/L. The period is 1 flux quantum for all 
states and the amplitudes are = 0.207, = 0.245 and TJ*^ = 0.268. 



0.3 -, 




-0.3-1 



exactly 1. We do not see any anomalous oscillations in the k = 2, 3,4 PF states, 
such as fractional flux periodicity of the persistent currents [10,31,32] which are 
characteristic for the BCS paired condensates, or more generally for some broken 
symmetries for all numerically accessible temperatures 0.03 < T/Tq < 14. This is 
a confirmation of the Bloch-Byers-Yang theorem [22] formulated in the context of 
superconductors. The persistent currents for the PF states have been shown to have 
oscillations for M = [32], however, this case is physically irrelevant. The case 
M = has only been used in [33] as a technical tool for investigating the wave 
functions for the PF states. One of the objectives of this paper was to show that the 
periods of the persistent currents for the PF states are exactly 1 for all temperature. 
Nevertheless, for T = it is possible that the QH systems undergo phase transitions 
to some BCS-like phases where the periods could be less then 1 due to some broken 
symmetries. 
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Furthermore, the amplitudes of the persistent currents in the parafermion FQH 
states decay exponentially with temperature, as shown in Fig. 3. 

Fig. 3. Temperature decay of the persistent current's amplitudes in the k - 2,3 and 4 
parafermion FQH states (without the contribution from the two v = 1 Landau levels with 

opposite spins) computed numerically in units evp/L, for temperatures measured in units 
of Tq. The zero temperature amplitude in these units are Vk/2, i.e., 1/4,3/10 and 1/3 re- 
spectively. 




4.3 Low-temperature asymptotics of the persistent currents in the PF^t states 

Because the low-temperature asymptotics of the chiral partition function (44) in- 
cludes only universal parameters such asvn, Qqh and the proper quasihole energy 
Cqh the low-temperature behavior of the persistent current which we analyze in this 
section would have a universal form and would be valid for all FQH states. 

Taking the logarithm of Eq. (44) and differentiating with respect to (p like in Eq. (68) 
we get the following expression for the persistent current for \(p\ < 1/2 

I(T, cp) ^=^^ ^ i^-VH(f> + 2(2qhe-S^ sinh (27:^^ Q^h^jj • (69) 
4.3.1 The persistent current's amplitude for T = 

For r = the second term in the curly brackets in Eq. (69) vanishes which ex- 
presses the fact that for zero temperature the free energy is determined by the 
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ground state energy. Then the persistent current is simply proportional to the flux 
7(r = O,<^) = -VH^0, for -i<,^<i 

which corresponds to the so called saw-tooth curve. It riches its maximum at = 
-1/2 which is 

^max = T V// — for 7 = 0. 



We stress that the zero temperature behavior of the persistent currents is completely 
determined by the CZ factors (36). 



4.3.2 The persistent current's amplitude for <T <^Tq 

Because the persistent current is an odd periodic function of the flux we are going 
to consider the half -interval only -1/2 < < 0. Writing the sinh in Eq. (69) into 
its exponential form and ignoring QX^{2n^TQQ^(f)IT) which vanishes in the limit 
r — > we get 

liT,^) = ^{-v.<^-eqhe-S^e--^^^.^4. 

T^Tq Li \ ) 



The solution of dI(T, 0)/50 = for - 1 /2 < < gives the position of the current's 
maximum 



qh 



<Pmax(T) ^= -— — + 



1 



T«To 2n^Qqk To 



In 



^qh iQ 



Vh 



(70) 



and substituting Eq. (70) into Eq. (69) we get the following expression for the 
persistent current's amplitude 



r«ro L 



^qh 



1 



Qqh 2n^Qqh To 



1 +ln 



/J 



(71) 



Using the charge-statistics relation Eq. (46) and the expression for the proper quasi- 
hole energy Eq. (45) we can rewrite Eq. (71) in the form of Eq. (47) according to 
which the persistent current's amplitude decays logarithmically with increasing the 
temperature due to the thermal activation of quasiparticle-quasihole pairs. 



32 



4.4 High-temperature asymptotics of the persistent currents in the Zjt parafermion 
FQH states: extracting the universal decay exponents 

Unlike the analysis of the low-temperature asymptotics of the persistent currents in 
Sect. 4.3 which has a universal form depending only on the quasiparticle charges 
and proper quasihole energies the high-temperature asymptotics, which is also uni- 
versal, depends on the explicit form of the modular S matrix so in the next subsec- 
tions we shall illustrate it on the examples of the parafermion FQH states for 
k = 1,2) and 4. We show on Fig. 4 the logarithmic plot of the persistent currents 
amplitudes for the Fermi liquid with v = 1, the Laughlin state with y = 1/3 and 
the parafermion FQH states with k = 2,3, and 4 computed numerically from the 
corresponding CFT partition functions for temperature 0.03 < T/Tq < 14 which 
includes also the high-temperature regime T :^ Tq. 




10 -| 1 1 1 1 1 1 1 1 1 1 1 — ^ ' 1 

2 4 6 8 10 12 14 

Fig. 4. Logarithmic plot of the temperature decay of the persistent current's amplitudes for 
the Fermi liquid (vh - 1), the Laughlin state (vh - 1 /3) and the A: = 2, 3 and 4 parafermion 
FQH states 



4.5 k=2 



According to Eq. (57) the irreducible representations of the PF2 model are labelled 
by the pairs (/,p) where -1 < I <2 and p = 0, 1 with the restriction p > l-p mod 2. 
Using the explicit S matrix for the Pfafiian FQH state from Ref. [34] (see Eq. (5.8) 
there) we find the universal coefficients necessary for the computation of the high- 
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temperature limit of the partition function (50) 

' V2 + 1 for / = 0, p = 



i'=-i p'>i'-p' 



V2-1 for/ = 0, p=l 
otherwise 



For convenience the characters of the ^ = 2 parafermion FQH state are written 
explicitly in Eq. (A.l) and the chiral partition function (58) after S transformation 
takes the form 



Z2(T, 0) = ( V2 + i)xo.oir', n + i^- i);ro,i(T', n = 

- ( V2 + l) W, 2^ ; 8) + ( V2 - l) W, 2C; 8), 



(72) 



where we have dropped out a ^'-independent term proportional to the central charge, 
substituted ch(AQ + Ag)(T') - q'^, ch(A^ + Aj)(t') - q'^'^ and ignored q'^'^ as com- 
pared to q'^ since q' ^ when T ^ oo. Next, again dropping the r] function and 
taking only the leading terms in the K functions, i.e., n = for Ko(S) and n = 0, -1 
for ^4(8) we get 



Z2(T,(f>) - (V2+1) l-f2 



V2- r 
^+1, 



k'cos(20 



Substituting the T and <p according to Eq. (48) and Eq. (49), using the approxima- 
tion ln(l -I- X) ^ X for Z «c 1 and taking the derivative with respect to according 
to Eq. (68) we get the following high-temperature asymptotic expression for the 
persistent current (note that In^ksTQlcpQ = evpjL) 



h{T, <f>) - -ir'iT) sin(27r0), I^^iT) = ^ - 



(73) 



4.6 k = 3 



For this case Eq. (57) tells us that the irreducible representations are labelled by 
(l,p) where -2 < / < 2 and < p < 2 with the restriction p > I - p mod 3. The 
S matrix has been written explicitly in Eq. (5.7) in Ref. [8] and the sums of the 
S -matrix elements are 



^"=E E 



(i',P') 



l'=-2 p'>l'-p' 



-| (sin (f ) + 3 sin (f )) for / = 0, p = 
^(3sin(|)-sin(f))for/ = 0,p = 2 
otherwise 
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The 10 independent characters in this case are shown in Eq. (A.2) and the partition 
function (58) for = 3 after S transformation becomes 

Z,(T, ct>) = Vo(t', n + F''^Xo,2(r', O - 

- Koir', 3r ; 15) + F^'^q'^l') + 

+ Ks(t', 3^; 15) (f''\'" + F''Y"') + 
+ K_s(t', 3r ; 15) (fO'V'^' + F'^'Y'"') , 

where we have used ch(A„ + Ao)(t') ^ q'^\ ch(Ao + A^)(t') = ch(A^ + A^Xt') ^ 
q'^l^\ ch(Aj + Aj)(r') = ch(A2 + A2)(t') ^ q'^'^ and ch(Aj + A^W) ^ q'^'^ (using 
Eq. (56) for the CFT dimensions), again skipping the central charge term in the 
characters (60). Keeping only the lowest powers in q' 0, ignoring the Dedekind 
function and taking only n = in the ^ functions we get 



To get the persistent current in the high-temperature limit we apply the same scheme 
as for ^ = 2 which yields 



h{T,<l>) - -7^^(r)sin(27r0), 

T»To 



5 2_ 
3 ^ 15 



(74) 



where we have used that 
/r0,2 _ 1 



/ro,o £, + r 



where D = 2 cos 



g) 



is the quantum dimension of the elementary quasihole in the PF3 FQH state [8]. 



4.7 k = 4 



For k = 4 the representations are labelled by (see Eq. (57)) (l,p) where -2 < / < 3 
and < p < 3 with the restriction p> l-p mod 4. The 10 dimensional S matrix S 
for the (neutral) affine coset is shown in Eq. (B.2) and the complete 15 dimensional 
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S matrix is computed by Eq. (65). Then the sums of the 5 -matrix elements are 



l'=-2 p'>l'-p' 



2 + V3 for / = 0, p = 

2 - V3 for 1 = 0, p = 2 

1 for 1 = 0, p = 3 
otherwise 



so that the partition function (58) for ^ = 3 becomes 



+ Ke(T', 4r ; 24) + F^'^q'^'"'^'''' + F^'^q'^'" + 

+ K_e(T', AC; 24) (f^'^^'^'^^'^^ + f''^'^""^'''^ + F^'^q""""^'''^) , 



where the complete characters are shown in Eq. (A.3) and we have used ch(A^ + 

A^)(t') ^ q'^'^^'P\ with APP(0,0) = 0, A^^(0,2) = 1/12, A^^(l, 1) = A^^(3,3) = 
3/4, APP(1,3) = 1/3, APP(2,2) = 1 (see Eq. (56) for the CFT dimensions A^^), 
again skipping the central charge term in the characters (60). Keeping only the 
lowest powers in q' 0, ignoring the Dedekind function and taking only n = in 
the K functions we get 

Z,(T,cf>) - F0'0 + 2F0'3^'^"^"^°''^cos(27rr) 

T»To 



To get the persistent current in the high-temperature limit we apply the same scheme 
as before 



h(T,(f>) ^-^^ -C^(r)sin(27r0), 



3 2_ 
2 12 



(75) 



Finally, we can summarize the results for the high-temperature asymptotics of the 
persistent current's amplitude for the PF^t FQH states 



rmax 



(76) 
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where the exponent ak for the Zjt parafermion FQH states is universal and can be 
written in general as 



<.. = i-+2A"(A„ + Aj. 2A"(A„+A,) = ?|^. (77) 

Remark 4 While a = \ jvufor the Fermi and Laughlin states [ 3,35 ], this exponent 
is not determined by the filling factor alone in the parafermion FQH states. 
Instead, there is a crucial contribution from the neutral sector which is described 
here by the parafermion coset. 

The analysis of the high-temperature asymptotics allows us to predict the form of 
the universal exponents a for arbitrary FQH states. Because F^ only if the 
corresponding electric charge Q{X) vanishes the neutral contribution to a would be 
determined from the neutral weight with minimal neutral CFT dimension which 
satisfies the pairing rule ^ with / = du, i.e., 

or = — + 2A^")(AJ, where P[AJ = dn mod A: with A("^(AJ = min,(78) 



i.e., A^ is the lowest weight of the neutral character having the minimal CFT dimen- 
sion among all neutral characters which couple to the m(1) character ^±d^(T, n^^; nHdn). 
Thus Eq. (77) should be valid for the entire parafermion hierarchy since A^ = 
A,) + A2 where A^^(A() + A^) = is minimal dimension of the weights satis- 
fying the pairing rule /'[A^ + A2] = 2 = k + 2 mod k. Accidentally, the parafermion 
CFT dimension A^^(Aq + A2) = for ^ = 2 (see Eq. (56)) and the exponent in that 
case is given by the inverse of the filling factor only. 

Using MAPLE we have computed the values of a and /q for the k = 2, 3, 4 PF states 
from the numerical data shown on Fig. 4 by linear regression (after subtracting the 
ln(r/ro) contribution) over 80 points in the interval 6 < T/Tq < 14. The exact 
values of the exponents a and the (interpolated) asymptotic amplitudes Iq for the 
first several PF^ states are in perfect agreement with the numerical ones and both 
are summarized in Table 1 . The exponents a/, for the parafermion FQH states have 
to be compared with those corresponding to Fermi and Laughlin liquids which 
according to Remark 4 are or = 1 and or = 3 respectively. As mentioned before, 
Eq. (76) demonstrates that at high temperature there exists a diff'erent mechanism 
reducing the persistent currents - thermal decoherence which leads to exponential 
decay of the persistent currents amplitudes. 



Note that the characters K^^j^ always appear in the vacuum character xo,o together with 
Kq and have the minimal nonzero CFT dimension 
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Table 1 

The slopes a and the j-intercepts Iq (in units evf/L) of the logarithmic plots in the high- 
temperature region for the A: = 2, 3, 4 PF states. 



k 


a 


/o [evf /L] 


numeric 


exact 


numeric 


exact 


2 


2.000706 


2 


0.110253 


(2/;r)(V2-l)/(V2+l) 


3 


1.800010 


1.8 


0.150309 


(2/7r)(2cos(;r/5) - l)/(2cos(;r/5) + 1) 


4 


1.666306 


5/3 


0.169827 


(2/;r)/(2 + V3) 



5 Conclusion 



We have shown that the two-dimensional CFT for the edge states of a FQH disk 
sample can be effectively used to compute analytically the partition function in 
presence of Aharonov-Bohm flux. The effect of the AB flux is interpreted as twist- 
ing of the electron operator and of the CFT Hamiltonian leading to the very com- 
pact formula (34) for the partition function in presence of flux. Using this partition 
function as a thermodynamical potential one can compute the oscillating persis- 
tent current, AB magnetization and magnetic susceptibility in arbitrary FQH states. 
Because of the absence of complete spin-charge separation the excitations of the 
FQH liquid satisfy a pairing rule (22) that we have formulated in a general 
form. As a result we find that in general the CFT characters should have a specific 
Znjj invariant structure (23) of sum of products of charged and neutral characters. 
The general approach has been illustrated on the example of the parafermion FQH 
states in the second Landau level using the results of Ref . [8] for the corresponding 
CFT. Our analysis shows that the persistent currents in the parafermion FQH states 
are odd periodic functions of the AB flux with period exactly one flux quantum as 
it should be according to the Bloch-Byers-Yang theorem. We find that these cur- 
rents are thermally suppressed by two diff'erent mechanisms in the low- and high- 
temperature limits, respectively and have universal non-Fermi liquid behavior. 

For the high-temperature asymptotics of the parafermion CFT we have computed 
explicitly the modular S matrices (65) using the convenient quasiparticle basis of 
Ref. [8]. As a byproduct we have derived the S matrix (66) for the diagonal affine 
coset (53) in a very compact form and have proven that it is equivalent to the S 
matrix of the su{k)2 current algebra, i.e., not only the irreducible representations of 
the coset (53) and su(k)2 are labelled by the same weights but the fusion rules of 
both CFTs are identical. 

The universal exponents (77) extracted from the high-temperature asymptotics of 
the persistent currents show that for the FQH states with numerator of the filling 
factor hh > I there is always neutral contribution depending on the neutral sector 
and the l^nH pairing rule. The experimental measurement of these exponents may be 
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useful for the identification of the correct CFT describing a given FQH universaUty 
class. 
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A Chiral partition functions for the parafermion QH states 

For ^ = 2 the 6 independent characters entering Eq. (58) using the notation Ki (8) := 
Ki (r, 2^; 8) and ch(ju,p) := ch (a^ + A^j (t) with < ;u < p < 1 can be written 

X-2fl{T, = K-2 (8) ch (0, 0) + K2 (8) ch (1, 1) 
;r~2,i(T, = K-2 (8) ch (1,1) + K2 (8) ch (0, 0) 
X-i,i(T, = K-1 (8) ch (0, 1) + K, (8) ch (0, 1) 

Xofiir, = Ko (8) ch (0, 0) + (8) ch (1, 1) 

Xo,i(t, = Ko (8) ch(l,l) + K, (8) ch (0, 0) 

XuiiT, = K, (8) ch (0, 1) + K_^ (8) ch (0, 1) (A. 1) 

For the ^ = 3 parafermion state we use the following notation: /T/ (15) := Ki (r, 3^; 15) 
and ch(//,p) := ch (A + A ) (r) with < < p < 2 the 10 inequivalent characters 
are 

X-2,i = K_2 (15) ch (0, 1) + (15) ch (1, 2) + K_j (15) ch (0, 2) 
X-2,2 = K_2 (15) ch (2, 2) + K^ (15) ch (0, 0) + K_, (15) ch (1, 1) 
X-u = (15) ch (0, 2) + (15) ch (0, 1) + (15) ch (1, 2) 
X-i,2 = K_, (15) ch (1,1) + K, (15) ch (2, 2) + K_e (15) ch (0, 0) 
ATo.o = ii^o (15) ch (0, 0) + iSTs (15) ch (1, 1) + is:_5 (15) ch (2, 2) 
;ro,2 = i^o (15) ch (1 , 2) + ^5 (15) ch (0, 2) + (15) ch (0, 1) 
Xu = K, (15) ch (0, 1) + (15) ch (1, 2) + K , (15) ch (0, 2) 
;ru = i^i(15)ch(2,2) + i^6(15)ch(0,0) + i^_4(15)ch(l,l) 
X2,i =if2(15)ch(0,2) + i^7(15)ch(0, 1) + if_3(15)ch(l,2) 
Xi,2 = K2 (15) ch (1,1) + (15) ch (2, 2) + K.^ (15) ch (0, 0) (A.2) 
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For k = 4 parafermion state where Ki (24) : = Ki (r, 4^ ; 24) and ch(jx, p) : = ch (a^ + 4p) W 
with 0<//<p<3we have 15 independent characters 



X-3,1 = (24) ch (0, 1) + K, (24) ch (1, 2) + K9 (24) ch (2, 3) + (24) ch (0, 3) 

X-3,3 = K 3 (24) ch (2, 3) + K3 (24) ch (0, 3) + K, (24) ch (0, 1) + K., (24) ch (1 , 2) 
X-2,i = K 2 (24) ch (0, 2) + i^:4 (24) ch (1, 3) + K^, (24) ch (0, 2) + (24) ch (1, 3) 
X-2,2 = K.2 (24) ch (1 , 1) + is:4 (24) ch (2, 2) + K,,, (24) ch (3,3) + i^_8 (24) ch (0, 0) 
X-2,3 = K 2 (24) ch (3,3) + ^4 (24) ch (0, 0) + (24) ch (1,1) + (24) ch (2, 2) 
X-i,2 = K., (24) ch (0, 3) + K5 (24) ch (0, 1) + (24) ch (1,2) + (24) ch (2, 3) 
X-i,3 = (24) ch (1,2) + (24) ch (2, 3) + Kn (24) ch (0, 3) + K_j (24) ch (0, 1) 
Xo,o = Ko (24) ch (0, 0) + Ke (24) ch (1, 1) + Kn (24) ch (2, 2) + K_e (24) ch (3, 3) 
Xo,2 = Ko (24) ch (1,3) + Ke (24) ch (0, 2) + /sTij (24) ch (1, 3) + K_e (24) ch (0, 2) 
Xo,3 = Ko (24) ch (2, 2) + isTe (24) ch (3, 3) + Kn (24) ch (0, 0) + K_e (24) ch (1, 1) 
Xu = K, (24) ch (0, 1) + Kj (24) ch (1, 2) + isT.n (24) ch (2, 3) + K_, (24) ch (0, 3) 
Xi,3 = K, (24) ch (2, 3) + is:7 (24) ch (0, 3) + (24) ch (0, 1) + K ^ (24) ch (1, 2) 
X2,i = K2 (24) ch (0, 2) + Ks (24) ch (1, 3) + ^_i„ (24) ch (0, 2) + K , (24) ch (1, 3) 
X2,2 = K2 (24) ch (1,1) + K, (24) ch (2, 2) + i^.io (24) ch (3, 3) + K ^ (24) ch (0, 0) 
X2,3 = K2 (24) ch (3, 3) + K^ (24) ch (0, 0) + K_,o (24) ch (1,1) + K_, (24) ch (2, 2) 

(A.3) 



B The diagonal coset's S matrix 



In this appendix we shall sketch the derivation of the modular S matrix for the coset 
(53) following the general approach of Goddard-Kent-Olive [26]. The characters 
of the diagonal coset appear by definition as string functions in the decomposition 
of the coset numerator su(k)i ® su(k)i in terms of the su(k)2 characters 

/l\r)/:^(r) = J] P(A, A'; A") ch(^^,f )(T);^g(T), (B.l) 

A" ^ ' 



where and;^f^^ are the su(k)i characters at level / = 1 and 2 respectively and the 
function 



P(A,A';A") = 



1 if [A"] = [A] + [A'] 
otherwise 



expresses the Zjt charge conservation rule [17]. Let us now apply the S transfor- 
mation r ^ -l/r to both sides, i.e., SxT = Z.S^iy^K S/l^ = Z^^'^iV^ 
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"^^A" = Zo-'^A''o--^^?^ Schp" = 'Zp"Sp»chp" and then decompose again the 
product ;(fj/^(T);f^!\T) in the left hand side in tenns of the;^^^^,^ and the coset charac- 
ters using Eq. (B.l) 



mod k fi" ' ' 

= J] P(A. A'; A") ^ ( "'l: ° ) W-f^W' 

A" p" a ^ ' 



where we have used the triple equivalence to label the coset representations by the 
level-2 index only. Next we substitute cr = fj." in the RHS, compare the coefficients 
of;f^^,\ multiply both sides by (S^^}^)~^ and sum over ju" to get 



' mod k fi" ^ ' 

P(A, A'; v)2^;"^h(^;;!; 



def 

Because (jx, ; fi") = (A^, A^, ; A^, + Ag) we can shift the summation variables in 
the LHS according to 

fi=Ji+IJ.', a = a + iJ.', = + fi' 



and use the triples equivalence (A^, A^,; A^ + Ag) = (A~, 0; A~ + Ag) to get for the 
LHS 

fx' mod k n mod k a</3 



where we have used that P(fl + fi' ,fi';a + fi' + fi + fi') = P(fi, 0; or + fi). Now we can 
use the properties of the S matrices under the action of simple currents 



c(l) _ p2;^^•^^^r.(l) ^(1) _ Iml^^H) ^(2) _ M r, (2) 



where Q^^^, (v) is the monodromy charge (62) in the su(k)2 CFT. Writing the su(k)2 
weight V = A^ + A^ with 0<p<o-<k-l and using the following formula for 

the CFT dimensions in su(k)2 

. (2). . , A N = 2p(/r -o-) + ik+l) [pjk -p) + o-jk - o-)] 
^-P 2k(k + 2) 
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we get the following simple expression for the monodromy charge 



e^i'CA +A.) = - 



k 



The final step is to take the summation over jj! using again the 5-function ^ 



set p" = A5 + (then jx = K^„^ ='• [p"] due to the conservation of the charge 

in the coset) and compare the coefficients of ch| ^ | ' ° j(T) which gives the coset S 
matrix, Eq. (66). 

For the high-temperature asymptotics of the persistent current in the ^ = 4 para- 
fermion FQH state we shall need the (neutral) S matrix (66). Using the following 
basis of 10 coset weights 

Ai + Ai + A2 A2 A3 + A3' 

the matrix S has been computed with Maple by Eq. (66) where the S matrix of 
the denominator su{A)2 of the coset has been computed by the Kac formula (67) 



IJ.' mod k 




V6V2V6V2V6V2V6V6V2V6 
12 4 6 4 12 4 6 12 4 12 



^ i±£ — -1+t A - V2 -1-t -/V2 

44^4 4 4^44 4 




44^4 4 4^44 4 




_^ /V2 We -i V2 -V6 -i V2 Ve V6 rv^ -V6 
12 4 6 4 12 4 6 12 4 12 



(B.2) 



-1+' V2 l+i n -V2 1-i i'VI 

44^4 4 4^44 4 



6 6 6 6 6 6 



_^ W2 _^ V6 _^ V6 W2 V6 

12 4 6 4 12 4 6 12 4 12 



-1-i rj -l+i iV2 1-i rv -V2 l+i -iVI 
44^4 4 4^44 4 



-iV2 -V6 iVI -V6 iVI V6 V6 -iV2 -V6 
12 4 6 4 12 4 6 12 4 12 



^ note that the A:-ality of v is [v] = [A + A^] - p + a 
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